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We report on a first principles theory for analyzing tlie 
parametric electron pump at a finite frequency. The pump 
is controlled by two pumping parameters with phase differ- 
ence (f). In the zero frequency limit, our theory predicts the 
well known result that the pumped current is proportional to 
sin (j). For the more general situation of a finite frequency, our 
theory predicts a non-vanishing pumped current even when 
the two driving forces are in phase, in agreement with the 
recent experimental results. We present the physical mecha- 
nism behind the nonzero pumped current at = 0, which we 
found to be due to photon-assisted processes. 
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Parametric electron pumpOlII is an interesting device 
which deUvers a finite DC current to the outside world at 
zero bias potential, by cyclic variations of two device con- 
trol parameters. Recently, an adiabatic quantum electron 
pump was reported in an open quantum dot where the 
pumping signal was produced in resfionse to the cyclic de- 
formation of the confining potentiaH. It was found that 
the pumping signal, Vdot (</>), is sinusoidal in the phase 
difference (j) between the two deforming potentials in the 
weak pumping regime, and it becomes non-sinusoidal in 
the strong pumping regime. The amplitude of the pump- 
ing signal increases linearly with the frequency of-the de- 
formation. Most notable, however, was the datap show- 
ing Vdot(O) ^ significantly at strong pumping, whereas 
ydot{T^) ~ for all pumping strength. Even in the weaker 
pumping regima. small deviations from Vdot(O) = could 
already be seenB. Thej-.tijaditional and successful para- 
metric pumping theoryB'cl, valid in the adiabatic regime 
and up to first order in frequency, requires two pumping 
parameters which traverse in a closed path in parameter 
space in each cyclic period and the pumping signal is pro- 
portional to the enclosed area by the path. Accordingly, 
if the two pumping parameters are in phase so that the 
enclosed area is zero, the pumping signal Vdot(O) should 
vanish. It is extremely puzzling that the experimental 
dataa consistently showed Vdot(O) ^ 0, a fact which has 
not been understood so far, and, clearly, it calls for the 
development of a first principle theory going beyond the 
adiabatic regime and low frequencies. 

It is the purpose of this letter to present a theory for 
parametric pumping which is valid at finite frequency. 
Using this theory we investigate the frequency depen- 
dence of the pumping current Ip{4>) — G dotVdoti4>) with 
Gdot a constant, and it allows us to understand why it is 
possible to have a pumping signal even when the pump- 
ing forces are exactly in-phase. When frequency is low. 



our theory recovers the traditional resultsElu. As the fre- 
quency increases, we predict a nonzero /p(0) which is a 
consequence of photon assisted processes, and our the- 
ory also predicts /p(7r) sa 0. These results allows us to 
reach the conclusion that the experimental dataa show- 
ing Vdot(O) 7^ and Vdoti'T'') ~ 0, are generic nonlin- 
ear transport features of parametric pumps at finite fre- 
quency. Furthermore, our theory predicts that at very 
large frequencies /p(7r) should start to deviate from zero; 
and most intriguing is the natural theoretical outcome 
that only one periodic deforming potential can produce 
a pumping signal at finite frequency. The latter is due to 
the fact that a finite frequency provides extra degrees of 
freedom through photon assisted processes that is capa- 
ble of playing the role of a second pumping parameter. 
These new predictions should be testable experimentally. 
We start by considering a parametric pump which con- 
sists of a coherence quantum scattering region attached 
to two ideal leads L, R. The leads maintain identical elec- 
trochemical potentiaL|?.e., fiL = fJ-R = M- The Hamilto- 
nian of this system isu 
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+ {Ud+d,+i + c.c.)] + Y,{Tk,cC+^d, + c.c.) (1) 



where Cka and di are annihilation operators of elec- 
trons for the jo. lead and the scattering region at site 
i, respectivelya. The three terms describe the leads, the 
scattering region, and the the coupling between the leads 
and the scattering region with the hopping matrix T^.q. 
In the last term, j = 1/A^ for a = L/ R. The parametric 
pump works by cyclic deformations of potential at two 
different pumping sites i and j in the scattering region, 
Vi/j{t) = Vi/j cos{ujt + (pi/j) with ipi/j the phase of the 
pumping force. 

Neglecting interaction between electrons in the ideal 
leads, the Keldysh nonequilibrium Green's function 
(NEGF) theory gives the following standard expression 
for the time dependent currentE£l {h — 1), 



(2) 
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where the Green's functions G'"'"'^ and the self-energy 
Sq are defined as in the usual mannertil. 

It is 

tedious but straightforward to evaluate tift-Green's func- 
tions by iterating the equation of motiont30, and for our 



purpose it is adequate to calculate them to second order 
in the pumping potential l^/j (t) which gives the average 
pumped current to the second order. We note that the 
next higher non-vanishing order is the fourth. After all 
the Green's functions arc obtained, the average pumped 
current is calculated from Eq. (^ by integrating time over 
one pumping cycle. We obtain, 
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where 



iqV^Vj 



ygr^G^rcona-/-) 



X e'^G"!;- - G^;-] + (/ - /+)e-^*[G^;+ - G^;+]} (4) 

where / = f{E) and f± = f{E ± w) are the Fermi dis- 
tribution functions; Tl — — 2/m[S£] is the line width 
function; (f> = ipj — ipi is the phase difference between 
the two pumping forces. Finally, G"'' = G^^{E) and 
QGr± ^ G'^^{E ± w) are the retarded Green's functions 
when there is no pumping force. In Eq.ffl), la is ob- 
tained by setting j — i. Eqs.(^,||) are the main result of 
this work. 

Before we discuss Eqs.(||,^) in connection to the ex- 
perimental data of Ref. ||, let's first examine the low fre- 
quency limit of these results and show that the_nonven- 
tional parametric pumping theory is recoveredd'u. We 
expand Eq.M) in powers of w and onlji_keep the linear 
term, this gives the adiabatic currentaQ. Note that la 
and Ijj are of the second order in frequency, Eq. (m re- 
duces to 
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whei£ we have used the relation G^ — G"" = —iG^TG"' 
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Using the Fisher-Lee 

reiation:JO S = -I + iT^^^GT^/^, Eq.(|) reduces to 
exactly the same exgression as that obtained from scat- 
tering matrix theorytl. Due to the factor sin(0) in the low, 
frequency result (0), one obtains the familiar outcomecl 
that Ip{(l)) = at both = and cj) = tt. 

Clearly, the low frequjency result (|^) does not explain 
the experimental resultH of /d(0 = 0) ^ 0, we need to in- 
vestigate the full result Eq. (Q) at a finite frequency. The 
first term on the right hand side of (^) has a clear physical 
meaning: it represents the photon absorption process in- 
dicated by the factor exp{iujt -\- iff)) . Similarly, the second 
term corresponds to the photon emission process with 
factor exp(— iwt — ifj)). These two competing processes 
tend to cancel each other in the expression of pumped 
current. These photon assisted processes are essential to 
understand the experimental finding that the pumped 
signal is nonzero at (f) — 0. To sec this clearly, we rewrite 
Eq.(0) into the following form 
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X [\gliE,E-cj) + e-'^^g]iE,E~u;)\' 

+ \gl'{E, E-Lo) + e-^^gf{E, E c.)p] 

+ U - U)[\9]{E,E + uo) + e^^g]{E,E + uj)\^ 

+ \g^{E,E + Lo)+e^^gf{E,E + u)\^]} (6) 



where gf{E,E~ lo) ^ ^r^TrZ[G\l{E)]V.^t{E - u:)] 
with a — 1,N indicating the left and right positions 
where scattering region is contacted by the leads, so that 
Fi = Fi and F^v = Tr- The propagator G5^„ describes 
the free motion of an electron from position n to position 
m in the device. Therefore, the quantity gl{E,E — uj) 
describes the following process: a charge carrier with en- 
ergy E~u! enters the device from the left lead, it absorbs 
a photon with frequency u at site i, and then exits from 
the left lead with energy E. This process is represented 
by the Feynman diagram in the inset of Fig.(l). Sim- 
ilarly, gj{E,E — uj) describes exactly the same process 
except that the electron absorbs a photon at position j. 
Now the physics is transparent: the first term in Eq.(@) 
represents an interference of the photon-absorption pro- 
cesses happening at positions i,j, i.e. the interference 
of the two processes in the Feynman diagram of Fig.(l). 
Therefore, when = 0, we have a constructive interfer- 
ence so that the g}{E, E—u) and the g]{E, E~uj) process 
add up (the factor exp(— i0) = -1-1). When = tt, there 
is a destructive interference in which the two photon- 
absorption processes cancel to a large extent (the fac- 
tor exp(— i0) = —1). Exactly the same can be said 
for the three other terms of Eq.(^): the term involving 
gf{E,E — oj) describes photon-absorption process with 
electrons entering the device at the right lead and exits 
from the left; and the terms involving gl{E, E + uj) and 
gf^{E, E + oj) describe photon-emission processes. 

In addition to the interference effects, the pumped cur- 
rent is also affected by a competition between photon 
emission and absorption, marked by the (/ — /_) term 
for absorption and the (/ — /+) term for emission, in 
Eq.(y). The combined effect, competition plus interfer- 
ence, can be clearly seen by expanding Eq.(p|) to order 
0{u>'^) which produces a complex expression for general 
0, but if we concentrate only on = or tt, the result is 
much simpler and physically transparent. 
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Ip = Y^^'i^E, - dE,) [\9UEi,E2) + e''^g]{EuE2 
-f \g^{E^,E2) + e^^gf{EuE2)\'] 
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where one sets Ei = E2 = E after taking the derivative 
and the temperature is set to zero. Hence, at = the 
constructive interference gives larger current than at = 
IT (see also the right inset of Fig. 2 described below). The 
behavior of the pumped current for general is the result 
of the interplay between the photon assisted processes 
and interference processes. 



The above physical picture allows us to unambiguously 
conclude that the behavior of pumped current depends on 
a combination of interference and competition of photon 
assisted processes. The competition between absorption 
and emission sets an overall magnitude for the pumped 
current at each 0. Destructive interference occurs at 
4> = ir, therefore Ipijr) « 0; constructive interference 
occurs at (j) = 0, giving rise to a non- vanishing lp{0). 
To demonstrate this, in the right inset of Fig. 2, we plot 
the magnitude of pumped current versus at a finite 
frequency due to the contributions of the photon absorp- 
tion (solid line, first two terms of Eq.(||)) and photon 
emission (dotted hne, the last two terms of Eq.(||)) pro- 
cesses, respectively. Taking into account of absorption 
and QMission processes, /p(0) is significantly larger than 
Ip (7r)o. The result in this inset is obtained from an one- 
dimensional model which we now discuss in more detail. 

Our discussion so far is completely general on the 
general result Eq.(Q). In the following we investigate 
a model in a more specific manner by applying Eq.(H) 
to a double barrier quantum structure which we model 
by potential U{x) = VoS{x) + VoS{x — a) whem Vq is 
the barrier height and a the barrier separationl3. For 
this system thft-jGreen's function G{x,x') can be cal- 
culated exactlyta. We choose the pumping force as 
V{t) = ViS{x-Xi)sm{u}t) + VjS{x-Xj)sm{ujt + (j)). With 
this specific pump modeL Eq. (0) can be evaluated numer- 
ically without difScultyU. 

In Fig.l, we plot the pumped current Ip at zero tem- 
perature versus the Fermi energy at different frequen- 
cies with phase difference (j) — tt/2. For comparison, we 
also plot transmission coefficient (long dashed line) . The 
peak in transmission coefficient indicates a quantum res- 
onance mediated by the resonance state in the double 
barrier. Clearly, the pumped current Ip also shows a 
resonance behavior, peaked at the same resonance state 
and is largely suppressed away from it. As the frequency 
is increased, the pumped current reverses the sign and 
the peak is shifted slightly. The general feature of Ip 
versus phase difference (jj can be obtained from Eq.(|j) 
which can be rewritten as Ip = ci + ci sin -I- C3 cos ~ 
c\ -I- C4 sin((?!) + 0o), where qs and (/)o are constants. This 
indicates that Ip has a sinusoidal behavior in (/). The 
sinusoidal behavior is a direct consequence of the fact 
that our theory is to the second order of the pumping 
amplitudetS. 

Fig. 2 presents Ip as a function of for several different 
frequencies at resonance where the transmission coeffi- 
cient is unity. We notice that at very low frequency (solid 
line) Ip ~ sin((/)). As the frequency increases, the ampli- 
tude of Ip also increases. At the same time, the entire 
/p-(/) curve is shifted upward although leaving /p((/' = tt) 
largely unchanged. However, /^(^ = 0) is seen to increase 
significantly for the curve with w = 0.13. As already dis- 
cussed above, these features are due to interference of 
the photon assisted processes, and they are ii^-very good 
agreement with the experimental observationa. The left 
inset of Fig. 2 shows Ip as a function of frequency for 



several phase differences. We observe that at small fre- 
quencies /p « at = and tt. For </) — -k jl and 
3/27r, Ip have similar values. At larger frequencies, Ip 
increases in different fashion. For instance, /p at (/) = 
increases slowly at small frequency and then has a lin- 
ear behavior for larger frequencies. At = 7r/2, Ip in- 
creases initially, reaches a maximum, and then decreases 
slowly. At (/) = TT, /p is very small but nonzero and in- 
creases linearly with a much smaller slope. Whereas Ip 
a,t (f) = 37r/2 is the largest among all the curves. This 
suggests that one should operate the parametric electron 
pump at (j) — 3tt/2 to achieve the maximum pumped 
current. 

So far, we have explained the physics behind the pe- 
culiar experimental finding of Ip{4> = 0) 7^ 0: it is due 
to photon assisted processes which is a nonlinear phe- 
nomenon. However, when the two pumping forces act 
exactly in-phase, nothing distinguishes them (except that 
they act at different positions of the pump). Therefore, it 
is extremely interesting to ask: can one operate a pump 
with only one pumping parameter ? To check if this is 
possible, we set V,- = in Eq.(3|) and notice that the 
resulting Ip still remains finite due to the first term of 
Eq.(||). Fig. 3 plots this Ip versus energy at different fre- 
quencies. Again, Ip is peaked near the resonant point. As 
the frequency increases, the amplitude of Jp also increases 
and its peak shifts towards larger energy due to photon 
assisted process. The inset of Fig. 3 shows Ip vs the po- 
sition of the single pumping site at different frequencies 
for an energy the resonance. Due to the symmetry of 
the system, Ip is antisymmetric across the center of the 
double barrier. It is surprising that the magnitude of 
the pumped current has the same order of magnitude as 
that pumped by two driving forces (see long dashed line 
in the inset of Fig. 2). The reason that a pump can oper- 
ate with only one external driving force is directly related 
to photon assisted processes which happen at a finite fre- 
quency. It will be interesting to test this prediction by 
further experimental work. 

In summary, we have developed a parametric pump- 
ing theory valid at finite frequency and it shows that 
the pumped signal for two pumping parameters in-phase 
can deviate from zero due to photon assisted process at 
finite frequency. This explains the anomaly at (/) = 
observed experimentally. This theory also suggests that 
even with one pumping parameter, it is possible to pro- 
duce a pumped signal at finite frequency whose ampli- 
tude is in the same order of magnitude of that of two 
pumping parameters. 

After the paper was submitted, we were made aware 
of a recent work by BrouweiO which examined magnetic 
field symmetry of the pump. It suggested that rectifi- 
cation effect might be important in understanding the 
= anomaly. 
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FIG. 2. /p versus (j) for different ui at resonant point. 
Solid line: uj = 0.01; dotted line: u> = 0.05; dash-dotted line: 
u! = 0.09; dashed line: uj — 0.13. Left inset: Ip versus ui 
for different values of (p at resonant point. Solid line: (j) = 0; 
dotted line: (j> = 7r/2; dash-dotted line: (p — n; dashed line: 
— 3tt/2. The long dashed line is Ip vs uj for single pumping 
parameter at position x — 0.2a. Right inset: the magnitude 
of pumped current at resonant point due to contributions of 
the first two terms (solid line) and due to the last two terms 
(dotted line) in Eq.(7), for uj = 0.13. 



FIG. 3. Ip versus energy at oj = 0.2, 0.4, 0.6, 0.8 near 
resonant point with a single pumping parameter. The po- 
sition of the pumping parameter is at a; = 0.2a. Inset: Ip 
versus position (in unit of a) of the pumping parameter at 
w = 0.2, 0.4, 0.6, 0.8. 



FIG. 1. Pumped current (solid line) and transmission co- 
efficient T (long dashed line) versus energy with (j> = 7r/2 at 
UJ — 0.2, 0.4, 0.6, 0.8. For illustrating purpose, we shifted T 
by multipling a factor 0.1. Inset: Feynman diagrams corre- 
sponds to Eq.(7). In Fig.l and Fig.2, we have set Xi = 0.02a, 
79.2, V^ = VJ = 1. 
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